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BACKGROUND OF THE INVENTION 



Field of the Invention 

[0001] The present invention relates generally to systems and methods for determining 
characteristics of a formation surrounding a borehole. More particularly, the present invention 
relates to an acoustic logging tool that provides enhanced performance through use of a quadrapole 
quadrupol e source configuration. 

Description of Related Art ^^"^^ 
[0002] Acoustic v^ell logging is a well-developed art, and details of acojjsrtclogging tools and 
techniques are set forth in A. Kurkjian, et al., "Slovmess^^Bgfimation from Sonic Logging 
Waveforms", Geoexploration, Vol. 277, pp. 215-2^6<r991); C. F. Morris et al., "A New Sonic 
Array Tool for Full Waveform Loggins,>^E-13285, Society of Petroleum Engineers (1984); A. 
R. Harrison et al., "Acquisition^,af{a Analysis of Sonic Waveforms From a Borehole Monopole and 
Dipole Source . . . "SPE 20557, pp. 267-282 (September 1990); and C. V. Kimball and T. L. 
Marzetta, "S^jaWance Processing of Borehole Acoustic Array Data", Geophysics, Vol. 49, pp. 
274^^*r(March 1984), all of which are hereby incorporated by reference herein, 
[0003] An acoustic logging tool typically includes an acoustic source (transmitter), and a set of 



receivers that are spaced several inches or feet apart. An acoustic signal is transmitted by the 
acoustic source and received at the receivers of the borehole tool which are spaced apart from the 
acoustic source. Measurements are repeated every few inches as the tool passes along the borehole. 
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[0004] The acoustic signal from source travels through the formation adjacent the borehole to the 
receiver array, and the arrival times and perhaps other characteristics of the receiver responses are 
recorded. Typically, compressional wave (P-wave), shear v^ave (S-wave), and Stoneley wave 
arrivals and waveforms are detected by the receivers and are processed. The processing of the data 
is often accomplished uphole or may be processed real time in the tool itself. Regardless, the 
information that is recorded is typically used to find formation characteristics such as formation 
slowness (the inverse of acoustic speed), from which pore pressure, porosity, and other formation 
property determinations can be made. In some tools, the acoustic signals may even be used to 
image the formation, 

[0005] Logging-while-drilling (LWD) logging tools are generally located as close to the drill bit as 
possible, so as to minimize the delay between reaching a formation and measuring its properties. 
When implemented as LWD tools, acoustic logging tools must overcome a number of obstacles to 
perform successfiiUy. These obstacles include drilling noise, and acoustic properties of the thick 
tool body. The drilling operation creates a continuous low-frequency noise that can interfere with 
acoustic measurements. For the most part, the noise of the drilling operation exists only in a range 
of frequencies below about 2 kHz. Accordingly, the drilling noise can be effectively screened out 
in tools designed to operate in higher frequency ranges. 

[0006] The portion of the drill string near the drill bit is typically designed to withstand the large 
compressive force created by the weight of the drill string resting on the drill bit. Accordingly, the 
walls of the tubing in this region is substantially thicker than the average wall thickness of the drill 
string. The outer-diameter of the tubing is consequently increased, particularly for LWD tools, 
which have to accommodate the bulk of logging instruments in addition to an inner bore for fluid 
flow. 
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MAMMID) W STOOTTUTIE SPIECIMCATIION 

[0007] In acoustic tools, increasing the rigidity of the tool increases the amplitude of "the tool 
mode", i.e. the amplitude of the acoustic energy that propagates through the body of the tool 
between the acoustic source and the receiver array. It is desirable to minimize the tool mode 
because this energy can interfere with the desired measurements of the true formation shear wave 
velocity. 

[0008] Increasing the outer diameter of the tool (relative to the borehole diameter) also increases 
the dispersion of interface waves. Interface waves are acoustic energy that propagates along the 
borehole boundary. For soft formations where the shear wave velocity is slower than the borehole 
fluid sound velocity, the interface wave velocity provides the best measurable indication of the true 
formation velocity. Dispersion of these waves reduces the measurement accuracy. Accordingly, it 
is desirable to reduce the dispersion of interface waves. 

SUMMARY OF THE INVENTION 
[0009] Accordingly, the rgjs disclosed herein an acoustic tool that provides a reduced tool mode 
and enhanced accuracy for est imating she ar wave propagation slowness in slow formations. In one 

: « ■ ._ 

embodiment, the acoustic tool comprises; an acoustic source, an array of acoustic receivers, and an 
intemal controller. The acoustic source excites waves that propagate in a quadrupole quadrupol e 
mode. The intemal controller processes signals from the array of acoustic receivers to determine a 
peak phase semblance having a slowness value that varies with frequency. The minimum slowness 
value associated with the peak phase semblance provides an accurate estimate of the shear wave 
propagation slowness. The acoustic source preferably includes four source elements. The elements 
that are 90° apart are preferably driven in inverse-phase to obtain the quadrapol e quadrupol e 
excitation pattern. 
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MAMIKEB-ILJF SILJIBOTTOTIE OTCIDMCATjIOM 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0010] For a more detailed description of the preferred embodiment of the present invention, 
reference will now be made to the accompanying drawings, wherein: ^ 

Fig. 1 shows a typical logging while drilling environment; 

Fig. 2 shows a preferred acoustic logging tool embodiment; 

Fig. 3 shows a cross-sectional view of the acoustic source portion of the logging tool; 
Figs. 4A-4C show monopole, dipole, and quadrapole quadrupol e excitation pattems; 
Figs. 5 A-5D show simulation results for monopole excitation; 
Figs. 6A-6D show simulation results for dipole excitation; 
Figs. 7A-7D show simulation results for quadrapol e quadrupol e excitation; and 
Fig. 8 shows a flowchart of a method for determining shear wave propagation slowness 
in slow formations. 

[0011] While the invention is susceptible to various modifications and alternative forms, specific 
embodiments thereof are showoi by way of example in the drawings, and will herein be described in 
detail. It should be understood, however, that the drawings and detailed description thereto are not 
intended to limit the invention to the particular fomi disclosed, but on the contrary, the intention is 
to cover all modifications, equivalents and alternatives falling within the spirit and scope of the 
present invention as defined by the appended claims. 

[0012] It is noted that speed and slovraess are inversely related, and the determination of one 
automatically determines the other. The following discussion focuses on determining wave 
propagation slowness, but this is recognized to be the equivalent of determining wave propagation 
speed. 
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MAMSEB-=I[JF STOmTHJTIE MCmCATON 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
[0013] Turning now to the figures, Fig. 1 illustrates a typical logging-while-drilling environment. 

* 

A drilling platform 2 is equipped with a derrick 4 that supports a hoist 6. The drilling is carried out 
by a string of drill pipes connected together by "tool" joints 7 so as to form a drill string 8. The 
hoist 6 suspends a kelly 10 that is used to lower the drill string 8 through rotary table 12. 
Connected to the lower end of the drill string 8 is a drill bit 14. The bit 14 is rotated and drilling 
accomplished by rotating the drill string 8, by use of a downhole motor near the drill bit, or by both 
methods. Drilling fluid, termed "mud", is pumped by mud recirculation equipment 16 through 
supply pipe 18, through drilling kelly 10, and down through the drill string 8 at high pressures and 
volumes to emerge through nozzles or jets in the drill bit 14. The mud then travels back up the hole 
via the annulus formed between the exterior of the drill string 8 and the borehole wall 20, through a 
blowout preventer (not specifically shown), and into a mud pit 24 on the surface. On the surface, 
the drilling mud is cleaned and then recirculated by recirculation equipment 16. The drilling mud is 
used to cool the drill bit 14, to carry cuttings fi-om the base of the bore to the surface, and to 
balance the hydrostatic pressure in the rock formations. 

[0014] In a preferred embodiment, downhole sensors 26 are coupled to a downhole 
controller/telemetry transmitter 28 that transmits telemetry signals by modulating the mud flow in 
drill string 8. A telemetry receiver may be coupled to the kelly 10 to receive transmitted telemetry 
signals. Other telemetry transmission techniques are well known and may be used. One of the 
sensors 26 is an LWD acoustic tool located near the drill bit 14. 

[0015] Fig. 2 shows a preferred embodiment of an acoustic LWD tool 40 in a borehole 20. The 
tool comprises an acoustic source 42, an acoustic isolator 44, and an acoustic receiver array 46. 
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The acoustic isolator 44 serves to attenuate and delay acoustic waves that propagate through the 
body of the tool from the source 42 to the receiver array 46. Any standard acoustic isolator may be 
used. The acoustic source 42 and receiver array 46 are preferably made up of piezoelectric 
elements designed to operate in downhole conditions. However, many other sources and detectors 
are suitable for downhole operation, and may be used. 

[0016] The LWD tool 40 preferably includes an internal controller (not specifically shown) that 
communicates with the downhole controller 28. When the downhole controller 28 enables the 
operation of the LWD tool 40, the intemal controller controls the triggering and timing of the 
acoustic source 42, and records and processes the signals from the receiver array 46. The intemal 
controller fires the acoustic source 42 periodically, producing acoustic pressure waves that 
propagate through the fluid in borehole 20 and into the surrounding formation. At the borehole 
boundary, some of the acoustic energy is converted into shear waves that propagate along the 
interface between the borehole fluid and the formation. As these "interface waves" propagate past 
the receiver array 46, they cause pressure variations that can be detected by the receiver array 
elements. The receiver array signals are processed by the intemal controller to determine the tme 
formation shear velocity. The intemal controller may then communicate the formation shear 
velocity to the downhole controller/telemetry transmitter 28 for communication to the surface. 
Altematively, the downhole controller may maintain a log of formation shear velocities dovrahole. 
[0017] Fig. 3 shows a cross-sectional view of the LWD tool 40. The acoustic source 42 preferably 
consists of four elements that are equally spaced on the circumference of the tool. The receiver 
array preferably includes a row of elements spaced axially along the surface of the tool. The 
receiver array may also include additional rows of elements spaced circumferentially from the first 
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row. In a preferred embodiment, the receiver array consists of four rows of elements that are 
circumferentially aligned with the four source elements. 

[0018] The source elements can typically be fired individually and concurrently. When fired 
concurrently, the elements can be fired in phase and in inverse phase. An in-phase firing means 
that the same trigger signal is applied simultaneously to the elements. An inverse-phase firing 
means that as one element receives a trigger signal, another receives an inverted trigger signal, i.e. 
the additive inverse (negative) of the trigger signal. This firing flexibility allows for variation of the 
acoustic source excitation pattem. 

[0019] Fig. 4A shows the four source elements being fired in-phase. This produces a monopole 
excitation pattem. This monopole pattem is one mode of wave propagation. Other modes of wave 
propagation are shown in Fig. 4B (dipole) and Fig. 4C ( quadrapole quadrupol e). The source 
produces a dipole excitation pattem if elements on opposite sides of the tool are fired in inverse- 
phase. A quadrapol e quadmpol e excitation pattem is produced when elements 90° apart are fired in 
inverse phase. 

[0020] Figs. 5-7 show simulation results for these three source excitation patterns. The simulations 
were performed using a concentric symmetric layer model. In these simulations, the formation is 
assumed to be a typical slow formation, with a compressional wave propagation slowness of 110 
/xs/ft, and a shear wave propagation slowness of 230 /is/ft. The borehole diameter is assumed to be 
8.5 inches. The inner diameter of the tool is 1.9 inches, and the outer diameter is 7.25 inches. The 
receiver elements are located at half-foot intervals between 4.0 ft and 7.5 ft, inclusive. Full 
waveform simulation has been done using existing techniques (such as 3D-finite difference 
numerical modeling). 
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[0021] Fig. 5B shows the simulated signals received by the receiver elements. Each of the 
waveforms is shown for a corresponding receiver as a function of time since the transmitter firing. 
(Note the increased time delay before the acoustic wave reaches the increasingly distant receivers.) 
After recording the waveforms, the logging tool typically normalizes the waveform so that they 
have the same signal energy. In Fig. 5B, each waveform is labeled with a percentage to indicate the 
original maximum peak amplitude of the signal relative to the maximum peak amplitude of the 
first signal. 

[0022] Fig. 5B also shows a graduated series of sloping lines to indicate the relative waveform 
delays to be expected for given slowness values. Slower waves (those having larger slowness 
values) take longer to reach the more distant receiver elements, and accordingly, their effect on the 
recorded waveforms is increasingly delayed for larger distances. 

[0023] To identify waves and their slowness values, the acoustic logging tool may calculate the 
time semblance E{t^) as a function of slowness and time for the data. This information in turn may 
be used to determine various formation properties, including wave propagation velocity and 
dispersion of acoustic waves. The equation for the time semblance E{tyS) is: 



1 



(N Y 



(1) 



In the above equation, N is the number of receiver elements, and hence is also the number of 
recorded waveforms, Xi{t) is the waveform recorded by the ith receiver, di is the distance of the zth 
receiver from the transmitter, and s is the slowness, hi Equation 1, the quantity (t-sdi) is the relative 
time at the zth receiver for a given slowness 5. 
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[0024] Semblance values E{t^) range between zero and one. Values near one indicate a high 
correlation between the various recorded waveforms at the given time and slowness, and hence 
indicate the presence of a propagating wave having that slowness value. Values near zero indicate 
little correlation between the various waveforms at the given time and slowness value, and hence 
provide no indication of a propagating wave having that slowness value. 

[0025] Fig. 5A shows the time semblance E(t,s) plot for the data in Fig. 5B. Two strong peaks are 
visible, along with three secondary peaks. The strong peak in the lower left comer corresponds to a 
compressional wave with the correct slowness (110 /xs/ft). The strong peak in the middle of the 
plot, and the secondary peaks below and to the right of the strong peak, correspond to the shear 
wave. The correct slowness value (230 /is/ft) is provided by one of the secondary peaks, but the 
multiplicity and relative strength of the peaks makes it difficult to correctly determine the shear 
velocity. 

[0026] Another function that may be useful for analyzing the performance of acoustic logging 
tools is the phase semblance. The phase semblance E(f^) is determined as a function of slowness 
and frequency according to the following equation: 



1 



2 



(2) 



(=1 



In the above equation, is the number of receivers, and hence is also the number of recorded 
waveforms, Xi(/) is the Fourier transform of waveform jc,<0, di is the distance of the zth receiver 
from the transmitter, and exponential factor e"^^^^^^*^' is the Fourier transform equivalent of the 
relative time shift (t-sdi). <S> represents the phase operator of the complex number. For a complex 
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number Ae!^, the phase given by the phase operator is <S>{Ael^]=d, (No attempt is made to hmit the 
phase to a range of principal values.) 

[0027] Fig. 5C shows the phase semblance plot for the data in Fig. 5B. The phase semblance is 
shown as a function of frequency (in kHz) and slowness (in /is/fl). A general phase semblance 
relationship is visible between frequency and slowness (sloping downward from left to right), but 
no clear indication of the shear wave is visible. 

[0028] Fig. 5D shows a plot of a term in the analytic calculation of the dispersion function as a 
function of frequency and slowness. The model used for this calculation is axi-symmetric. 
Standard continuous stress-displacement boundary conditions are used, along with a vanishing 
shear stress in the fluid. The dispersion function is found by predicting the behavior of the different 
wave propagation modes as a function of frequency. The techniques for performing these 
calculations may be found in the following references, which are hereby incorporated herein by 
reference: 

K.M. White, Underground Sound, Elsevier, Amsterdam, 1983 

D.P. Schmitt and M. Bouchon, Full wave acoustic logging: synthetic microseismograms 
and fi*equency wavenumber analysis, Geophysics, 50, 1756-1778, 1985, 
[0029] At any given frequency, four peaks are visible. The slowness value (and width) of the 
peaks varies as a function of frequency. The peaks are indicative of branch cuts and poles in the 
dispersion function. In general, the effect of the branch cuts on the waveform is negligible 
compared to that of poles. Two of the peaks represent branch cuts. These are relatively constant at 
1 GO /is/fl and 200 /is/fl, respectively. The bottom-most peak represents wave propagation along the 
interface between the tool and the borehole fluid, and the top-most peak represents wave 
propagation along the interface between the borehole fluid and the formation (i.e. the borehole 
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boundary). Most significantly, the propagation slowness of the borehole boundary wave gradually 
decreases as a function of frequency from 350 /xs/ft to about 280 /xs/fl. Because the true shear wave 
propagation slowness is 230 /xs/ft, there is little hope for an accurate measurement of the formation 
shear wave slowness using the monopole source configuration. 

[0030] Fig. 6B shows the signals received by the receiver elements for the dipole excitation 
pattem. These signals are the basis for the semblance plot in Fig. 6A and the phase semblance plot 
in Fig. 6C. The semblance plot in Fig. 6A shows one strong peak at about 800 lis and 150 /xs/ft, 
and three secondary peaks at about (1320 jits, 150 pts/fl), (1900 /xs, 250 /xs/ft), and (2400 /xs, 150 
/xs/ft), respectively. None of these peaks provides the correct slowness for the formation shear 
wave. 

[0031] The phase semblance plot in Fig. 6C shows two clear phase relationships which slope 
downward from left to right. The lower one is indicative of the tool interface wave, and the upper 
one is indicative of the borehole interface wave. The presence of both waves in the region of 
interest may cause one of the interface waves to be obscured by the other, which may complicates 
measurement procedures. 

[0032] Fig. 6D shows the analytic term for the dipole excitation pattem. Below 10 kHz, four peaks 
are apparent. (Additional peaks emerge at higher frequencies, but these are ignored for the present 
discussion.) As before, two of the peaks have relatively constant slowness values at 100 /xs/fl and 
200 /xs/ft, respectively. The tool interface wave starts at about 230 /xs/fl and falls off asymptotically 
to about 100 /xs/ft. The borehole interface wave varies between 300 /xs/fl and 280 /xs/fl. As before, 
there is little hope for obtaining an accurate measurement of the formation shear wave slowness 
from the borehole interface wave. 
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MAMKIEIID^IIJF SOTOTTOTIE MCfflMCATIdDN 

[0033] Fig. 7B shows the signals received by the receiver elements for the quadrapole quadrupol e 
excitation pattern. These signals are the basis for the semblance plot in Fig. 7A and the phase 
semblance plot in Fig. 7C. The semblance plot shows a single strong peak at 2000 (is and 275 /xs/ft, 
no help is apparent here. The phase semblance plot, however, shows a clear phase relationship 
between the slowness and the frequency. Significantly, the phase semblance peak has a sharp 
minimum at 240 /xs/ft. Using the quadrapol e quadrupol e excitation pattem and the phase 
semblance plot provides an improved method for determining the formation shear slowness. 
[0034] Fig. 7D provides additional support for the use of the quadrapol e quadrupol e excitation 
pattem. The analytic term for calculating dispersion has the famiUar four peaks. The peak 
associated with the borehole interface wave converges to the desired 230 /is/ft at about 2 kHz. For 
frequencies above 5 kHz, the peak maintains a slowness value of about 280 /xs/ft. 
[0035] From the above discussion, it is clear that the quadrapole quadrupol e source provides more 
direct results for determining the slow formation shear wave slowness than the monopole and 
dipole sources. A second substantial benefit of the quadrapol e quadrupol e source is that the tool 
interface wave is absent from the frequency region of interest. This eliminates one potential source 
of interference with the measurements, and simplifies the formation measurement procedure. 
[0036] Fig. 8 shows a flowchart of the improved method for determining formation slowness, hi 
block 52, the internal controller triggers the acoustic source in a quadrapole quadrupol e excitation 
pattem. hi block 54, the internal controller records the signals from the receiver array elements, hi 
block 56, the internal controller calculates the phase semblance as a function of frequency and 
slowness, and in block 58, the intemal controller identifies the minimum slowness value of the 
phase semblance peak. The intemal controller communicates this value to the downhole controller 
as the formation shear wave slowness. 
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[0037] Numerous variations and modifications will become apparent to those skilled in the art 
once the above disclosure is fully appreciated. It is intended that the following claims be 
interpreted to embrace all such variations and modifications. 

[0038] For example, the present method may be used with any acoustic source that excites wave 
propagation in the quadrapole quadrupol e mode, although a quadrapole quadrupol e source (i.e. a 
source that primarily excites the quadrapole quadrupol e mode) is preferred. The receiver 
configuration may be additionally or altematively altered (or the receiver signals processed) to 

I 

enhance the quadrapol e quadrupol e mode response and reduce the response to monopole and 
dipole modes. In one embodiment, the receiver array quadrapol e quadrupol e response may be 
enhanced by determining a combined signal for each set of four receivers in the array. The 
combined signal may be the sum of inverted signals fi-om two opposing receivers with the signals 
firom the remaining two receivers in the set. 
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